Identifying the neutrino mass hierarchy 
with supernova neutrinos 



M. KachelrieB and R. Tomas 



. M ax- Planck- Institut fiir Physik (Werner-Heisenberg-Institut), Miinchen 

(N 
o 

Q 



Abstract 



We review how a high-statistics observation of the neutrino signal from a future 
galactic core-collapse supernova (SN) may be used to discriminate between differ- 
ed ■ ent neutrino mixing scenarios. Most SN neutrinos are emitted in the accretion and 
cooling phase, during which the flavor-dependent differences of the emitted neutrino 
I spectra are small and rather uncertain. Therefore the discrimination between neu- 
■ trino mixing scenarios using these neutrinos should rely on observables independent 
. of the SN neutrino spectra. We discuss two complementary methods that allow 
CN ! for the positive identification of the mass hierarchy without knowledge of the emit- 
' ted neutrino fluxes, provided that the 13-mixing angle is large, sin^ i^is ^ 10~^. 
. These two approaches are the observation of modulations in the neutrino spectra 
by Earth matter effects or by the passage of shock waves through the SN envelope. 
Q^' If the value of the 13-mixing angle is unknown, using additionally the information 
• encoded in the prompt neutronization burst — a robust feature found in all modern 
D , SN simulations — can be sufficient to fix both the neutrino hierarchy and to decide 



whether ■i?i3 is "small" or "large." 



H ■ 1 Introduction 



Despite the enormous progress of neutrino physics in the last decade, many open ques- 
tions remain to be solved. Among them are two, the mass hierarchy — normal versus 
inverted mass spectrum — and the value of the 13-mixing angle "i^ia, where the observation 
of neutrinos from a core-collapse supernova (SN) could provide important clues [HEIIS]- 
Schematically, the neutrino emission by a SN can be divided into four stages: Infall phase, 
neutronization burst, accretion, and Kelvin- Helmholtz cooling phase. Most SN neutrinos 
are emitted during the last two stages, in all flavors with only small differences between 
the Ue and u^^r spectra. Moreover, the absolute values of the average neutrino energies as 
well as the relative size of the luminosities during the accretion and cooling phases cannot 
be determined with sufficient precision, especially if the SN is optically obscured and the 
progenitor type remains unknown. Therefore, a straightforward extraction of oscillation 
parameters from the bulk of the SN neutrino signal seems hopeless. Only features in the 
detected neutrino spectra that are independent of unknown SN parameters should be used 
in such an analysis. 
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The two most promising sources for such features are the modulations in the neutrino 
spectra caused by the Earth matter or by the passage of shock waves through the SN 
envelope. In the first case, matter effects on SN neutrinos traversing the Earth give rise 
to specific frequencies in the energy spectrum of these neutrinos, which are analytically 
known and depend only on the neutrino properties and the distance traveled through the 
Earth [H IH] . In the other case, the passage of the SN shock waves through the density 
region corresponding to resonant neutrino oscillations with the atmospheric neutrino mass 
difference imprints specific time- and energy-dependent modulations on the neutrino energy 
spectrum |3 |Hj, difficult to be mimicked by other effects. Only the amplitude of both 
modulations, and thus the statistical confidence to detect them, depends on how different 
the emitted neutrino fluxes are, while the specific shape of the modulations is independent 
from the fluxes. 



Case 


Hierarchy 


sin^ -dis 


Earth 


Shock 


z/e burst 


A 


Normal 


> 10-' 


Yes 


No 


No 


B 


Inverted 


> 10-3 


No 


Yes 


Yes 


C 


Any 


< 10-^ 


Yes 


No 


Yes 



Table 1: The presence of Earth-matter and shock wave effects in the z/g spectra and of the 
z/e burst for different neutrino mixing scenarios. 

In the following sections we will concentrate on three different neutrino mixing schemes 
(A, B, C), cf. Tab. Q where modulations by Earth or SN shock effects are clearly 
separated. For an inverted hierarchy and intermediate values of the 13-mixing angle, 
10~^ < sin^-i^is < 10~', both effects can be present. In this case, the value of has to be 
close to 10~^, the lower limit of the intermediate range. In Fig. ^ we show schematically 
the sensitivity of several SN observables to different neutrino mixing parameters. 

If at the time of the SN detection the value of t^ia is known to be "large," then the 
neutrino mass hierarchy can be identified observing the modulations induced either by the 
SN shock wave propagation (case B in Tab. or by the Earth matter effects (case A). 
If the value of ^^13 is still unknown and Earth matter effects are observed, an ambiguity 
exists between case A and C. In sec. 4, we discuss why the neutrino emission during the 
neutronization burst is only weakly dependent on variations in the input of current SN 
models and how the neutronization burst can be used to break the degeneracy between A 
and C. 

2 Identifying signatures of the SN shock wave prop- 
agation 

The neutrino spectra F^^ arriving at the Earth are determined by the primary neutrino 
spectra as well as the neutrino mixing scenario, Fy^{E, t) = J2jPji{E, t)F^^{E, t), where 
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Figure 1: Sensitivity of tlie tliree SN observables discussed in the text to different neutrino 
mixing scenarios. 



Pji is tlie conversion probability of a Uj into z/j after propagation tlirougli tlie SN mantle. 
Tlie probabilities pji are basically determined by the number of resonances that the neu- 
trinos traverse and their adiabaticity. Both are directly connected to the neutrino mixing 
scheme. In contrast to the solar case, SN neutrinos must pass through two resonance layers: 
the H-resonance layer at pn ~ 10^ g/cm^ corresponding to Aml^^, and the L-resonance 
layer at pl ~ 10 g/cm^ corresponding to Am^. Whereas the L-resonance is always adi- 
abatic and in the neutrino channel, the adiabaticity of the H-resonance depends on the 
value of ^^13, and the resonance shows up in the neutrino or antineutrino channel for a 
normal or inverted mass hierarchy respectively |[T]. 

During approximately the first two seconds after core bounce, the neutrino survival 
probabilities are constant in time and in energy for all three cases A, B, and C. However, 
at t ^ 2 s the H-resonance layer is reached by the outgoing shock wave, see the left panel of 
Fig. 121 The way the shock wave passage affects the neutrino propagation strongly depends 
on the neutrino mixing scenario: cases A and C will not show any evidence of shock 
wave propagation in the observed z/g spectrum, either because there is no resonance in the 
antineutrino channel as in scenario A, or because the resonance is always strongly non- 
adiabatic as in scenario C. However, in scenario B, the sudden change in the density breaks 
the adiabaticity of the resonance, leading to observable consequences in the spectrum. 

The key ingredient to observe signatures of the shock wave propagation is the time and 
energy dependence of the neutrino survival probability. In the right panel of Fig.|21 we show 
p{E, t) = pjy^jy^ averaged with the energy resolution function of Super-Kamiokande, for the 
case with a forward and a reverse shock. The latter forms when a neutrino-driven baryonic 
wind develops and collides with the earlier, more slowly expanding SN ejecta. Although 
the exact propagation history depends on the detailed dynamics during the early stages of 
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Figure 2: Left: Shock and reverse-shock propagation. The density profile is shown at 
the indicated instances after core bounce. The density region pn corresponds to resonant 
neutrino oscillations with the atmospheric mass difference, pl to the solar one |H1. Right: 
Survival probability p{E,t) as function of energy at different times [S]. 

the SN explosion, a reverse shock forms in all models which were computed with sufficient 
resolution The presence of two shocks results in a dip in p[E,t) at those energies for 
which the resonance region is passed by both shock waves. All these structures move in 
time towards higher energies, as the shock waves reach regions with lower density. 

A useful observable to detect effects of the shock propagation is the average of the 
measured positron energies, (Eg), produced in inverse beta decays Ue + P ^ n + e"*". In 
Fig. El we show (E^) together with the one sigma errors expected for a megaton water 
Cherenkov detector and a SN in 10 kpc distance, with a time binning of 0.5 s: Both panels 
contains the case that no shock wave influences the neutrino propagation, the case of only a 
forward shock wave and of both forward and reverse shock wave. The left and right panels 
show two different models for neutrino fluxes: Gl assumes different average energies of the 
emitted neutrinos, (Eq^Ux)) / {Eo{i^e)) = 1-2, and similar fluxes, ^o{^e)/^o{^x) = 0.8, while 
G2 assumes identical energy spectra, (Eq^Ux)) / {Eo^Ue)) = 1, and ^o{i^e)/^o{^x) = 0.5. 

The effects of the shock wave propagation are clearly visible, independent of the as- 
sumptions about the initial neutrino spectra. Moreover, it is not only possible to detect the 
shock wave propagation in general, but also to identify the specific imprints of the forward 
and reverse shock versus the forward shock only case. The signature of the reverse shock 
is its double-dip structure compared to the one-dip of a forward shock only. To study the 
dependence of the double-dip structure on the value of ^13, we show {E^) as function of 
time for different 13-mixing angles in the left panel of Fig. |3] Even for as small values as 



4 




Figure 3: The average energy of up ne~^ events binned in time for a static density profile 
(magenta), a profile with only a forward shock (red) and with forward and reverse shock 
(blue). The error bars represent 1 a errors in any bin, from Ref. 0. 



tan^ "(913 = 5 X 10~^ the double-dip is still clearly visible, while for tan^ -i^is = 1 x 10~^ only 
a bump modulates the neutrino signal. 

In the right panel of Fig. |3J we show the number of events binned in energy intervals of 
10 MeV as function of time for the case of a reverse shock. We can observe clearly how the 
positions of the two dips change in each energy bin. It is remarkable that the double-dip 
feature allows one to trace the shock propagation: Given the neutrino mixing scheme, the 
neutrino energy fixes the resonance density. Therefore, the progress of the shock fronts 
can be read off from the position of the double-dip in the neutrino spectra of different 
energy. Thus, the observation of shock wave effects does not only identify case B (inverted 
hierarchy, large ^13), but gives also access to physics deep inside the SN. 



3 Earth-matter effects 

During the first two seconds after post-bounce, during which roughly half of all neutri- 
nos are emitted, the dependence of the probability to reach the Earth on the neutrino 
energy E is very weak. However, if neutrinos cross the Earth before reaching the detec- 
tor, the conversion probabilities may become energy-dependent and induce modulations in 
the neutrino energy spectrum. These modulations may be observed in the form of local 
peaks and valleys in the spectrum of the event rate ctFP plotted as a function of 1/E. 
These modulations arise in the antineutrino channel only in cases A and C. Therefore its 
observation would exclude case B. This distortion in the spectra could be measured by 
comparing the neutrino signal at two or more different detectors such that the neutrinos 
travel different distances through the Earth before reaching them jH |U] . However these 
Earth matter effects can be also identified in a single detector OEl- 
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Figure 4: Left: Time dependence of {Ee) for a profile with a forward and reverse shock for 
several values of ian^di^ as indicated. Right: Number of events binned per energy decade 
as function of time for forward and reverse shock, from Ref. [Hj. 




Figure 5: Left: Realistic power spectrum from a single simulation jH]. Right: Comparison 
of P95 as a function of nadir angle 77 using a "floating cut" as discussed in Ref. P for a 32 
kton scintillator (SC) and a megaton water Cherenkov (HK) detector. 
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The net flux at the detector may be written in the form 

7 

= sin^ ^912^0 + cos^ dr^Fl + AF° ^ A, sin\hyl2) , (1) 

i=l 

where y is the "inverse energy" parameter y = 12.5 MeV/E, AF^ = (Fg—F^) depends only 
on the primary neutrino spectra, whereas the Ai depend only on the mixing parameters 
and are independent of the primary spectra. 

The last term in Eq. is the Earth oscillation term that contains up to seven ana- 
lytically known frequencies ki in y, the coefficients AF^Ai being relatively slowly varying 
functions of y. The first two terms in Eq. are also slowly varying functions of y, and 
hence contain frequencies in y that are much smaller than the k^. The frequencies ki are 
completely independent of the primary neutrino spectra, and can be determined to a good 
accuracy from the knowledge of the solar oscillation parameters, the Earth matter density, 
and the position of the SN in the sky jB] . The latter can be determined with sufficient pre- 
cision even if the SN is optically obscured using the pointing capability of water Cherenkov 
neutrino detectors [TU] . 

The power spectrum of N detected neutrino events is 

N 2 

J2 e'''^' . (2) 

i=l 

In the absence of Earth effect modulations, G{k) has an average value of one for k > 40. 
The region k < 40 is dominated by the "0-peak" , cf. Fig. which is a manifestation 
of the low frequency terms in Eq. (P). Identifying Earth effects is equivalent to observing 
excess power in G{k) around the known frequencies ki. The area under the power spectrum 
between two fixed frequencies k^i^i and fcmax is on an average (fcmax — ^min)- In the absence 
of Earth effects, this area will have a distribution centered around this mean. The Earth 
effect peaks tend to increase this area. The confidence level of peak identification, pa, may 
then be defined as the fraction of the area of the background distribution that is less than 
the actual area measured. 

In the right panel of Fig. El we assume the model Gl for the neutrino fiuxes and compare 
the results obtained with a 32 kton scintillator detector and a megaton water Cherenkov 
detector. In the latter neutrinos travel more and more distance in the mantle the 

peak moves to higher k values, and due to the high k suppression, the efficiency of peak 
identification decreases. When the neutrinos start traversing the core, additional low k 
peaks are generated and the efficiency increases again. 

The identification of Earth matter effects excludes case B, and is thus complementary 
to the observation of shock wave effects. 



G(k) = — 



4 Neutronization z/g burst 

If the value of is unknown, a degeneracy exists between case A and C. Both scenarios 
predict the same Ve signature in a water Cherenkov detector, and therefore the previous 
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Figure 6: Left: Neutrino luminosities as function of time for different progenitor masses. 
Right: Number of events per time bin for different reactions in a megaton water Cherenkov 
detector for a SN at 10 kpc for cases A (dashed hues) and C (sohd hues) and for different 
progenitor masses. Statistical errors are also shown for the 15 Mq case. From Ref. [11]. 

two observables are not useful to disentangle them. In this case, the additional information 
encoded in the Ve neutrinos emitted during the neutronization burst can fix the range of 
"i^is as well as the neutrino mass hierarchy. 

The prompt neutronization burst takes place during the first ~ 25 ms after the core 
bounce with a typical full width half maximum of 5-7 ms and a peak luminosity of 3.3- 
3.5x10^^ erg s~^. The striking similarity of the neutrino emission characteristics despite 
the variability in the properties of the pre-coUapse cores is caused by a regulation mecha- 
nism between electron number fraction and target abundances for electron capture which 
establishes similar electron fractions in the inner core during collapse. This leads to a 
convergence of the structure of the central part of the collapsing cores and only small 
differences in the evolution of different progenitors until shock breakout. The small depen- 
dence of the neutronization burst on, e.g., the progenitor mass can be verified in Fig. IHl 
left panel (cf. also Refs. [II1II2]). 

Theoretically, the identification of the neutronization burst is cleanest with a detector 
using the charged-current absorption of z/g neutrinos, the most abundant flavor during 
the burst. Examples of such detectors are heavy water detectors like SNO or liquid argon 
detectors like ICARUS [13] . The simplest possible observable to identify the neutronization 
burst is the total number of Ve events within an arbitrary fixed period tmax after the on- 
set of the neutrino signal. However, the probability that the SN is obscured by dust is 
as high as ~75%. Without an estimate for the SN distance, the total number of events 
observed cannot be connected to the SN luminosity and is thus not a useful observable. 
Instead, the time structure of the detected neutrino signal should be used as signature for 
the neutronization burst [TTj . 
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Since the event number in current and proposed charged- current detectors is not high 
enough to allow for a detailed time analysis, we discuss only the case of a megaton water 
Cherenkov detector. A draw-back of this choice is that this detector type does not have 
a clean signature for the channel. Instead, one has to consider the elastic scattering 
on electrons, which is affected by several backgrounds like inverse beta decay or reactions 
on oxygen. In Ref. it was shown that this background can be substantially reduced 
by using angular and energy cuts, as well as Gadolinium to tag neutrons from inverse beta 
decays. The sample of elastic scattering events still contains the irreducible background of 
scattering on electrons of other neutrino flavors than z/g, but this contamination does not 
affect the possibility to disentangle the different neutrino scenarios. 

The time evolution of the signal depends strongly on the neutrino mixing scheme. In 
case A, the survival probability is close to zero, and therefore the peak structure observed 
in the initial Vf. luminosity is absent. On the contrary, in case C, 30% of the original 
remain as whereas 70% are converted into v^. Since the cross section of z/g on electrons 
is much larger than that of z/j., the signal is dominated by the contribution of Vg. These 
z/g's follow the time evolution of Lj,^, and thus lead to a clear peak in the signal. 

In the right panel of Fig. IHlwe show the expected neutrino signal from t = —5 to 18 ms 
for different progenitor masses, and for the mixing scenarios A and C. The peak structure 
can be clearly seen in case C, but not in case A ^11^. Including recent improvements of the 
electron capture rates or uncertainties in the nuclear equation of state has only little effect 
on the neutronization peak compared to the size of the statistical fluctuations. Therefore 
the observation of a peak in the first milliseconds of the neutrino signal would rule out 
the normal mass hierarchy with "large" "i^is (case A), breaking the degeneracy between 
scenario A and C observed in the channel. 

After the neutrino mixing scheme has been established, the robustness of the theoret- 
ically predicted event number of the neutronization burst makes a measurement of the 
distance to a SN located at 10 kpc feasible with a precision of about 5%. Since it is 
likely that a Galactic SN is optically obscured by interstellar dust and no other distance 
determination can be used, this new method relying only on neutrinos is very promising. 

5 Summary 

A reliable determination of neutrino parameters using SN neutrinos should be indepen- 
dent from the primary neutrino fluxes produced during the accretion and cooling phase 
of the SN. Earth-matter effects and the passage of SN shocks through the H-resonance 
both introduce unique modulations in the neutrino energy spectrum that allow one their 
identification without knowledge of the primary neutrino spectra. While the observation 
of Earth-matter effects in the Ue energy spectrum rules out case B, modulations in the Ue 
time spectrum identify case B. If the value of ^9i3 would be known to be large, then the 
neutrino mass hierarchy would be identified. Otherwise, the detection of the neutroniza- 
tion z/g peak — a robust feature of all modern SN simulations — can break the remaining 
degeneracy between A and C. 
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